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Jet energy loss, photon production and photon-hadron correlations are studied together at high 
transverse momentum in relativistic heavy-ion collisions at RHIC energies. The modification of 
hard jets traversing a hot and dense nuclear medium is evaluated by consistently taking into account 
induced gluon radiation and elastic collisions. The production of high transverse momentum photons 
in Au+Au collisions at RHIC is calculated by incorporating a complete set of photon-production 
channels. Comparison with experimental photon production and photon-hadron correlation data 
is performed, using a (3-|-l)-dimensional relativistic hydrodynamic description of the thermalized 
medium created in these collisions. Our results demonstrate that the interaction between the hard 
jets and the soft medium is important for the study of photon production and of photon-hadron 
correlation at RHIC. 
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I. INTRODUCTION 

Hard processes are regarded as good tools for the to- 
mographic study of relativistic heavy-ion collisions and 
of the quark-gluon plasma: The energetic probes are cre- 
ated in the early stage of the collisions and can there- 
fore access the space-time history of the transient hot 
and dense nuclear medium created in these collisions. 
Large transverse momentum partons have been identi- 
fied as especially useful because they interact directly 
with the nuclear medium, and will lose energy in the pro- 
cess. The partonic jets which emerge will fragment into 
the hadrons which are observed in the detectors. Ex- 
periments at the Relativistic Heavy Ion Collider (RHIC) 
at Brookhaven National Laboratory (BNL) have indeed 
shown that high-pT hadrons in central A-l-A collisions 
are significantly suppressed in comparison with those in 
binary-scaled p-|-p interactions [J, 13, Hi • Those results 
are commonly referred to as "jet-quenching" . 

A lot of effort has been devoted to describe the en- 
ergy loss experienced by the color charges inside the hot, 
strongly interacting matter. Gluon bremsstrahlung with 
the Landau-Pomeranchuk-Migdal (LPM) [3| effect has 
been proposed as the dominant mechanism for the jet- 
quenching process, and different theoretical formalisms 
have been elaborated to describe the hard parton energy 
loss H 1, 1, 3., 1, 10,, lil, M, Several of these 

bremsstrahlun g ca lculations were compared against each 
other [3, [itI. Tisl. [Toj . using a relativistic ideal (3-1-1)- 
dimensional hydrodynamical simulation of the strongly 
interacting medium [20| . This implementation improved 
upon earlier attempts relying on simple schematic mod- 
els by building on a detailed and realistic description of 
the time evolution. Another potentially important en- 
ergy loss mechanism for high-px color charges is pro- 
vided by scattering off thermal partons, in binary elas- 
tic collisions. Estimations of the relative magnitude of 
coUisional versus radiative energy loss have also been 
performed in several different approaches and scenarios 
[2l|, m, m, m, in, 113, m . in Ref. m , a study of 



radiative and coUisional ener gy l oss was carried out in a 
given single approach ( AMY |l3l [Ti. [T5j ) . Note that the 
coUisional energy loss was further improved upon in Ref. 

Complementary to hadronic species, photons repre- 
sent another set of promising observables with the po- 
tential of probing the quark-gluon plasma (see, for ex- 
ample, |31| and references therein). Since they carry no 
color charge, once produced, they will essentially escape 
from the collision zone without further interaction with 
the surrounding nuclear medium, reflecting directly local 
medium properties. However, the experimental obser- 
vation of produced photons in high energy nuclear col- 
lisions involves the entire evolution of the collision, to- 
gether with a variety of sources. Considering that the 
background contribution from meson decay may in prin- 
ciple be experimentally subtracted, theoretical studies 
need to include photons produced in the early collisions, 
those involving jet-medium interactions during jet propa- 
gation in the medium |32| , and those from fragmentation 
of the surviving jets after their escape from the medium. 
In Ref. [H, [SJ], the calculation of photon production 
from nuclear collisions at RHIC with those different pho- 
ton sources has been done by first employing a (1-1-1)- 
dimensional Bjorken evolution model, and then (2+1)- 
dimensional relativistic hydrodynamics. 

These calculations, and others, have given rise to the 
hope that the "tomography" of the hot and dense core 
of the nuclear medium created in relativistic heavy-ion 
collisions may be put on a firm quantitative footing. 

In what concerns hadrons, several quantities have been 
put forward as "tomographic variables" . One of those 
is the nuclear modification factor Raa which is inferred 
from the measurement of single-particle inclusive px 
spectra. It is defined as the ratio of the hadron yield in 
A+A collisions to that in binary-scaled p+p interactions. 
Since the single-particle spectrum is a convolution of the 
jet production cross section and jet fragmentation func- 
tions, the suppression of the produced hadrons at a fixed 
Pt involves jets with a wide range of initial transverse en- 
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ergies; for details see [l^,[2^. In retrospect, it is therefore 
not particularly surprising that a wide variety of energy 
loss conjectures that differ significantly in the predicted 
energy loss mechanism can reproduce the measured Raa 
in central Au+Au collisions at ^ = 200 AGeV at RHIC 
given the present experimental errors for central colli- 
sions. In that sense, the tomographic usefulness of Raa 
in central collisions has become rather questionable [ssj . 
There have been several suggestions to improve the situ- 
ation. One is to consider a more differential observable, 
and study Raa in non-central collisions as a function of 
the azimuth and px [36j . This effectively corresponds to 
studying an average over different paths of partons as 
they traverse the expanding medium, which is no longer 
rotationally invariant in the transverse plane. In Ref. 
[l6t , Raa was also studied as a function of pT in central 
and non-central collisions at mid and forward rapidities. 
Although this does not restrict the initial jets' energies or 
initial vertices, it allows for a different averaging over the 
medium than in central collisions and can thus provide 
us with additional information. 

Another road to improvement consists of considering 
many-body variables, e.g. to study the production of 
high-pT liadrons correlated with other high-pT particles. 
One motivation is that the correlation between the trig- 
ger particle and the measured hadrons will constrain 
the initial energy distribution of the parton that frag- 
ments into the hadrons more tightly than single-particle 
measurements. One such suggested trigger is a hadron, 
see e.g., [13] ■ Note however that choosing a specific 
Pt for the trigger hadron on the "near-side" still does 
not strictly pin down the energy of the "away-side" par- 
tons. The trigger hadron is produced by fragmentation 
of partons that have also traversed the medium and lost 
energy. Therefore, a specific trigger bias is introduced 
into the away-side parton distribution and initial ver- 
tex distribution. This effect has been explored in [svj . 
where a detailed analysis of the initial vertex density and 
conditional probability distribution P{pt) of away-side 
partons for di-hadron correlations was presented in the 
Baier-Dokshitzer-Mueller-Peigne-Schiff (BDMPS) frame- 
work of energy loss Q by gluon brcmsstrahlung. There 
it was shown that the yield of hadrons per trigger hadron 
was sensitive to the specific assumptions about the evo- 
lution model for di-hadron correlations, while it was not 
for Raa- 

Another promising trigger is a high-p^ photon: One 
studies jet-quenching by measuring the pr distribution 
of charged hadrons in the opposite direction of a trig- 
ger direct photon (ssl . [39j . Direct photons are predom- 
inantly produced from hard collisions in the early stage 
of relativistic nuclear collisions. Assuming that a high- 
Pt direct photon from these processes can be used as a 
trigger of the away-side hadron, the transverse energy 
of the initial away-side parton before energy loss can 
then be deduced (up to corrections at next to leading 
order in a^). Thus, direct photons can provide a cali- 
brated probe for the study of the properties of high en- 



ergy density QCD. This kind of trigger does not intro- 
duce a vertex bias and therefore weighs all vertices ac- 
cording to the nuclear overlap function; and it strongly 
restricts the initial away-side partons' energy to a mono- 
energetic source. Given this restriction of the initial jet's 
transverse energy, the expected tomographic capabilities 
of photon-tagged jet correlation measurements arc much 
higher than those of Raa- It has been shown in Ref. [sHi 
that while different schematic conjectures of jet energy 
losses could not be discriminated by Raa in central col- 
lisions at mid-rapidity, they gave clearly different predic- 
tions for 7-hadron correlations. Some other studies have 
also been performed along this direction [iO, |4l|, \^ . 

However, jet-photon correlation calculations need to go 
beyond direct photons and to include the other important 
and known electromagnetic sources such as fragmenta- 
tion and jet-medium interaction. As the initial away-side 
partonic jets in these processes could have larger trans- 
verse energies than the near-side trigger photon, pho- 
tons produced from those processes might have signifi- 
cant contribution to the correlations between final pho- 
tons and hadrons. In this work, we consider these possi- 
bilities and incorporate a complete set of high-pr photon- 
production channels and study their relative contribu- 
tions to the correlations between back-to-back photons 
and hadrons. 

In Ref. [1^ , we have presented a consistent calculation 
of both coUisional and radiative energy loss in the same 
formalism. There it was applied to calculate the nuclear 
modification factor Raa of neutral pions in heavy ion col- 
lisions at RHIC for different centralities at mid-rapidity. 
Here, we continue and extend this effort by studying the 
effect of jet-medium interaction on the photon produc- 
tion as well as on correlated back-to-back hard photons 
and hadrons in high energy nuclear collisions. We cm- 
ploy the formalism developed in Ref. [2^ to account 
consistently for coUisional and radiative energy loss of 
the leading hard parton while it traverses the surround- 
ing soft nuclear matter, modeled by (3-|-f )-dimensional 
hydrodynamics [20| . 

This paper is organized as follows: the next two sec- 
tions describe our approach to jet energy loss and to 
photon production. Then, photon-hadron correlations 
are treated, and comparisons with experimental results 
obtained at RHIC by the PHENIX and STAR collabora- 
tions are performed. 



II. JET ENERGY LOSS 

In our approach, hard jets (quarks and gluons) evolve 
in the soft nuclear medium according to a set of Fokkcr- 
Planck type rate equations for their momentum distribu- 
tions P{E, t) = dN{E, t) /dE. The generic form of these 
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rate equations can be written as follows (33|, |43 



dt 



ab 



Pa{E + UJ,i) 



dTa^j{E + w, w) 



duodt 



dTj^b{E,Lj) 



diodt 



(1) 



where dTj^a{E,LL})/dLL}dt is the transition rate for the 
partonie proeess j — > a, with E the initial jet energy 
and Lo the energy lost in the proeess. The lo < part 
of the integration accounts for the contribution from the 
energy-gain channels. The radiative and collisional parts 
of the transition rates have been extensively discussed in 
Ref. [H, 111,13 • 

The initial jet momentum profiles may be computed 
from perturbative QCD calculations [Zsj . 



dcTAB- 



d^plpdy 



(2) 



ahd 



1 2XaXb daab^jd 
X : 

TT 2xa - xi^ev dt 



with = 2pip I y/SNN, where y/s^N is the center-of- 
mass energy. In the above equation, Ga/A{xa,Q) is the 
distribution function of parton a with momentum frac- 
tion Xa in the nucleus A at factorization scale Q, taken 
from CTEQ5 [461 including nuclear shadowing effects 
from EKS98 [43|. The distribution da/dt is the leading 
order QCD differential cross section with i^jot account- 
ing for next-to-leading order (NLO) corrections. It is 
fixed following Ref. [48| by reproducing the experimen- 
tal measurement of the inclusive 7r° cross section at high- 
Pt in in p-l-p collisions at ^/sNN = 200 GeV (see Fig. 1 
in Ref.[l^). The Kjot-factor is found to be 2.8 when 
the renormalization scale and the factorization scale are 
taken to be the transverse momentum of the initial jets, 
and the fragmentation scale is taken to be the transverse 
momentum of produced hadrons. 

To obtain the final high-px hadron spectra in Au-I-Au 
collisions at RHIC, the energy loss of partonie jets in 
the nuclear medium must be taken into account. This is 
performed by calculating the medium-modified fragmen- 
tation function D}^/j{z, r^, <j)) for a single partonie jet, 



Dh/3{z,r±, 



Dh^y{z')P{py\p,,r^,cl)),i3) 



where z = Ph/pj and z' ~ Ph/Pj' are two momentum 
fractions with ph the hadron momentum and pj {pj' ) the 
initial (final) jet momentum; Dfi/j{z) is the vacuum frag- 
mentation function, taken from the KKP parametriza- 
tion (49j . In the above equation, the probability function 
P{pj'\pj^ rj_, (f)) is obtained by solving Eq. ([T]), represent- 
ing the probability of obtaining a jet j' with momentum 
Pji from a given jet j with momentum pj. This depends 
on the path taken by the parton and the medium profile 
along that path, which in turn depends on the original lo- 
cation of the jet, rj_, and on its propagation angle with 



respect to the reaction plane. Therefore, one must con- 
volve the above expression over all transverse positions 
and directions to obtain the final hadron spectra: 



dcTAB- 



■hX 



1 



d^p!i,dy 27r 



d^fi_rABib,fi_) 



E 



d^pipdy 



dz ~ 

-i^,/,-(z,f^,0) — rr . (4) 



where VABih, r_\_) is the probability distribution of initial 
hard jets in the transverse position rj_, which is deter- 
mined from the overlap geometry between two nuclei in 
the transverse plane of the collision zone. For A-l-B col- 
lisions at impact parameter 6, 



rABib,r^_) 



TAir±+b/2)TBir^-b/2) 

TAB{b) 



(5) 



where TA(f'x) = J dzpA{r±, z) is the nuclear thickness 

function and TAB{b) = / d'^rj_TAirj_)TBir± + b) is the 
nuclear overlap function. A Woods-Saxon form for the 
nuclear density function p{r±, z) = po/[l -\- exp(^^^j^)] is 
used, and the values of the parameters R = 6.38 fm and 
d ^ 0.535 fm are taken from [soj . 

Putting all of the above together, one obtains the to- 
tal yield of hadrons produced in relativistic nuclear col- 
lisions. 



dN^Xsib) _ N,ou{b) d<JAB^ 



• hX 



d^plf^dy 



(JNN d^plldy 



(6) 



where TVcoU is the number of binary collisions and unn is 
the inelastic cross section of elementary nucleon-nucleon 
collisions. Finally, the nuclear modification factor Raa 
is defined as the ratio of the hadron yield in A-l-A colli- 
sions to that in p-l-p interactions scaled by the number 
of binary collisions. 



T)h 



AA 



1 dNl^{b)/d^pll,dy 
Ncoiiib) dNj^p/d^p'^dy 



(7) 



III. PHOTON PRODUCTION 

As mentioned already, there are several sources of high- 
Pt photons in relativistic nuclear collisions. The most 
important are direct photons, fragmentation photons, 
and jet-medium photons (bremsstrahlung photons and 
conversion photons). The thermal photon emission aris- 
ing from the partonie medium and later hadronic phase 
makes a negligible contribution in the high-pT regime, 
and thus is excluded from consideration in the present 
calculation. 

Direct photons are predominantly produced from early 
hard collisions between partons from the initial nuclei, 
via quark-anti-quark annihilation (q + q ^ g + "f) and 
quark-gluon Compton scattering {q{q) + g ^ q{q) + 7). 
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The inclusive direct photon cross section may be calcu- 
lated from Eq. ^ by applying the corresponding par- 
tonic processes for photon production, 



J^dircct 



E 



Q)Gb/B{xb,Q) 

abd 
1 2x^X5 



da. 



b — >-7tZ 



TT 2xa — xZey dt 



(8) 



with x7p 



2p^/y/sNN. The factor is a function of 
the photon's transverse momentum and is deduced by 
performing an NLO calculation of photon production in 
p-|-p collisions (sil . [s^ . [ssj , in which the renormalization 
scale and the factorization scale are set to be the photon 
transverse momentum. 

Fragmentation photons are produced by the surviving 
high energy jets after their passing through the hot and 
dense nuclear medium. The calculation of fragmentation 
photon spectra is similar to that of the high-py hadron 
production described in the last section. We may also 
define a medium-modified photon fragmentation function 
Dj/j{z, r±, <f>) as in Eq.Q. Then the final expression for 
the fragmentation photon cross section can be written as 



da 



frag 



d^pLdy 



^ j d^f^VAB{ri_) 



da A 



B^jX 



d^plpdy 



(9) 



where the vacuum fragmentation functions of quarks and 
gluons into real photons are taken from [53 |. 

Jet-medium photons are produced by jet-medium in- 
teractions during the passage of jets through the hot nu- 
clear medium. These include induced photon radiation 
(bremsstrahlung photons) and direct conversion from the 
high energy jets (conversion photons). It has been shown 
that those two processes are important for the under- 
standing of experimental data for pho ton production in 
Au+Au colhsions at RHIC [H,!!!,!!!. 

In order to incorporate the photons directly produced 
from jet-medium interaction, we may add to the set of 
evolution equations shown in Eq. ([1]) an additional evo- 
lution equation for photons. 



dP^^ 



dt 

where dT 



dujPqq{E+UJ,t) 



dTl%{E+u,u) 



diudt 



l^^/dujdt = dr^^^^^Vdwdt + dT'g°ji:^/dujdt. 



(10) 



The 



transition rates dV]^'^™ / dujdt for photon bremsstrahlung 
processes are similar to gluon bremsstrahlung processes 
and are extensively discussed in [H, [3, [11]. And the 
transition rates dT'^'^'^ / dujdt for binary jet-photon con- 
version processes may be inferred from the photon emis- 
sion rates for those processes [s^ , 



jpconv 



diodt 



E 



2 27ra,a,r2 



1 , ET 
— In — ^ 



3^; 



+ C2^2 (E/T) S{lo), (11) 



where = is the thermal quark mass and 

C2_2 « -0.36149 in the limit of S » T [H, [H, [13. 
The function S{lli) represents the fact that the produced 
photon has essentially the same energy as the incoming 
quarks (or anti-quarks) in the jet-photon conversion pro- 
cesses. 



IV. PHOTON-HADRON CORRELATIONS 

In this section, we present the calculation of the corre- 
lation between back-to-back hard photons and hadrons. 
The associated hadrons arc produced from the fragmen- 
tation of surviving jets after their passing through the 
nuclear medium, while the trigger photons may come 
from various sources as has been discussed in the pre- 
vious section. Therefore, the photon-hadron correlation 
will depend on the jet-photon correlation at the produc- 
tion vertex as well as on the energy loss of the jets during 
their propagation in the medium. 

In correlation studies, one of the commonly exploited 
observables is the yield per trigger, which is the momen- 
tum distribution of produced hadrons on the away side, 
given a trigger photon of momentum p^ in the near side. 
Following Ref. [58|, one may write 



Pipl) ' 



(12) 



where P{p^) = l/ctot j dy^da^ / dy^dplp represents 
the single-particle pr distribution and P{p^,pf^) = 
l/'^tot / dy^dyhda^+h/ dy^dyhdp'Jpdp!^ is the 7 - /i pair 
Pt distribution. 

The trigger photon and the associated hadron are pro- 
duced from a pair of back-to-back partons. Assuming no 
correlation between the individual evolutions of two par- 
tonic jets once they are produced, we may write down 
the expressions of these two distributions as follows 



Pf{pl) = I^Jd'f^VABir^) 



dp^P,{p=:,)P{pl\p^,r^, 



xE / dp'rpdplj,P^{p'j,,plj,) 



PfiPT.Pl) 



33' 

xP{pI\p^, fi, (/.)P(p^|y;, rl, 7r+<^). (13) 

In the above equations, Pf{plp) and Pf{p'^,p^) are the 
final state distribution functions for single-photon and 
back-to-back photon-hadron production, while Pi(plp) 

and Pi{pip,pip) are the initial distributions for a jet 
or jet pair to be produced in the medium with the 
given momenta and species types. The yield per trigger 
P{Pt\Pt^ ri, (/)) is the number of produced photons given 
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an initial jet, which may be decomposed into different 
contributions, 



0), (14) 



where the sources include direct photons, fragmentation 
photons and jet-medium photons as described in the last 
section. For the direct photon contribution, the near- 
side jet is replaced by a direct photon, so the yield 
per trigger is simply given by P(piJ,, directjpy, rx, 0) = 

^{Pt ~ Pt)- T^^^ fragmentation photon contribution is 
related to the medium-modified photon fragmentation 



function by P{p^, fragjp^, i 
with 



Dj/j{z,f±,(l))/pPj,, 

z = p7p/plp. Furthermore, we also need to in- 
clude the contribution from jet-medium photons, which is 
achieved by solving Eq. (fTOll . the photon evolution equa- 
tion, to obtain the yield per trigger P{pJp, JM|p^, r±, (j>). 
Note that for the case of hadron yield per jet trigger 
P{pt^\p'lp,f±,(j>), we only need to include the fragmenta- 
tion contribution since we focus on high-pj' hadron pro- 
duction. 

For the study of photon-hadron correlations, it is of- 
ten useful to define the photon-triggered fragmentation 
function as follows [H^, [60| , 

(15) 

with zt = p^/p^. In order to quantify the effect of 
the nuclear medium on the photon-hadron correlations, 
we may also introduce the nuclear modification factor 
Iaa^ which is defined as the ratio between A-l-A and p+p 
collisions of the photon-triggered fragmentation function. 



DAAizT^Pr) = PtPaa{Pt\Pt)^ 



Iaa{zt,Pt) 



Daa{zt,Pt) 
Dpp{zt,Pt) 



(16) 



V. RESULTS AT RHIC 

In the previous sections, we have outlined the theoret- 
ical formalism for calculating jet energy loss, the produc- 
tion of high-pT hadrons and photons, and photon-hadron 
correlations. In this section, we apply this methodol- 
ogy to study Au-I-Au collisions at RHIC, using a ( 3-1-1)- 
dimensional relativistic ideal hydrodynamics [20|, and 
compare our results to available experimental measure- 
ments. 

In Fig. [1] we show our earlier study [1^ of the neu- 
tral pion Raa in Au+Au collisions at RHIC measured at 
mid-rapidity for two different impact parameters, 2.4 fm 
and 7.5 fm, compared with the experimental measure- 
ments by PHENIX for the most central (0-5%) and 
mid-central (20-30%) collisions. In the calculation, the 
strong coupling constant as is adjusted in such a way 
that the experimental measurement of Raa in niost cen- 
tral collisions (upper panel) is described. The same value, 
as ~ 0.27, is used throughout the calculation. (There is 
no additional parameter for the later calculation of high- 
Pt photon production and photon-hadron correlations.) 
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FIG. 1: (Color online) The nuclear modification factor Raa 
for neutral pions in central and mid-central collisions (from 
Ref. [1^). The dashed curves account for only induced gluon 
radiation, the dash-dotted curves for only elastic collisions 
and the solid curves incorporate both elastic and inelastic 
energy losses. 



Fig. [T] also compares the relative contributions of induced 
gluon radiations and elastic collisions to the final Raa- 
One may find that while the shape does not show a strong 
sensitivity, the overall magnitude of Raa is sensitive to 
both radiative and coUisional energy loss. 
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FIG. 2: (Color online) Photon production in p+p collisions 
compared with PHENIX data [63,[6i|. 

Next, we present the results for high-pT photon pro- 
duction at RHIC. In Fig. [21 the inclusive photon cross 
section in p-f-p collisions at ^/snn = 200 GeV as a func- 
tion of photon Pt is compared with the experimental 
measurements by PHENIX jsil . [6^ . The theoretical cal- 
culation can nicely describe the experimental data; this 
serves as the baseline for calculating photon production 
in Au-|-Au collisions. Photons from early hard collisions 
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(Compton scatterings and annihilation processes) dom- 
inate at high-pT' regimes, while fragmentation photons 
gain increasing significance as photon pr decreases. Also 
shown is a power-law fit to the experimental measure- 
ment of total photon yield in p-l-p collisions at RHIC: 
da/dPpTdy = f .027/(1 + pr/0. 793 GeV)6**^3QgV-Vb. 
This function is somewhat closer to the low pT experi- 
mental data, and its usefulness will become clear shortly. 
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FIG. 3: (Color online) The contributions from different chan- 
nels to the photon production in Au+Au collisions at RHIC 
for h = 2.4 fm compared with PHENIX data [g^. 
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FIG. 4: (Color online) The contributions from different chan- 
nels to the photon production in Au+Au collisions at RHIC 
for b — 4.5 fm compared with date from the PHENIX collab- 
oration [gI, . 

In Fig. [21 the relative contributions from different chan- 
nels to high-pT photon production in central Au+Au 
collisions (6 = 2.4 fm) at RHIC is shown as a func- 
tion of photon pt and compared with most central 0- 
10% PHENIX data [H. While direct photons dominate 
the high-pT regime {pr > 6 GeV), the presence of jet- 
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FIG. 5: (Color online) The nuclear modification factor Raa 
calculated for photons in Au+Au collisions at RHIC for b — 
2.4 fm compared with PHENIX data with a 0-10% centrality 

cut fel. 



medium interaction is a significant contribution to the to- 
tal photon yield in Au+Au collisions at RHIC, especially 
in the intermediate-py regime {pt ~ 3-5 GeV). This can 
be more easily seen in Fig. 21 where we compare our 
calculation of photon production for 6 = 4.5 fm in the in- 
termediate Pt regime with 0-20% Au+Au collisions data 
from PHENIX At very low pr {pt < 2 GeV), 

the thermal emissions from partonic and hadronic phases 
are expected to dominate |64{ . This component is ex- 
cluded in the calculation as we are focusing on the high- 
Pt regime. Also the assumptions essential for jet-energy 
loss calculation break down at such low px- 

In order to further quantify nuclear medium effects on 
photon production in Au+Au collisions, we show in Fig. 
|5l the calculated photon Raa as a function of photon 
Pt for central Au+Au collisions {b = 2.4 fm) at RHIC 
compared with most 0-10% PHENIX data. The figure 
shows that Raa can be smaller than 1 even if only direct 
and fragmentation photons are included in both p+p and 
Au+Au calculations. As direct photons dominate the to- 
tal yield (and thus Raa) in the high-pj' regime, the de- 
crease of Raa at high pt is consistent with an isospin 
effect. If we include the contribtition from fragmentation 
photons in the calculation, photon Raa gets suppressed 
because these photons are produced from the fragmen- 
tation of the surviving jets with less energies due to the 
energy loss of jets when they arc traversing the thermal- 
izcd medium. However, the presence of jct-mcdium in- 
teraction in Au+Au collisions (but not in p+p collisions) 
again enhances Raa if we include them in the total pho- 
ton yield. It is telling that the sharp rise of Raa at 
low-pT originates mainly from the fact that our calcula- 
tion for p+p collisions is below the data points (see Fig. 
12]): a better agreement with the data is obtained with 
the power-law fit to the p+p data used to calculate Raa- 
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The additional QGP sources then manifest themselves 
in the difference between the solid and the dash-dotted 
curves in Fig. [5l The sensitivity of Raa to this level 
of details also makes clear the need for a precise, QCD- 
based, quantitative understanding of photon data in p+p 
collisions. 
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FIG. 6: (Color online) The photon momentum azimuthal 
asymmetry, V2, as a function of the photon transverse mo- 
mentum. The different components are explained in the text, 
and the data are from Ref. [6511. 
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FIG. 7: (Color online) The photon-triggered fragmentation 
function as a function of momentum fraction zt in p-fp col- 
lisions at RHIC. 

The nucleus- nucleus photon data can also be analyzed 
in terms of its anisotropy in momentum space, which 
is characterized by the well-known V2 coefficient. This 
data is show in Figure [H ft is clear that the high cen- 
trality cut will cause the momentum anisotropy to be 
small, and this is indeed what is observed. Moreover, the 
intrinsically negative V2 associated with jet-medium pho- 
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FIG. 8: (Color online) The photon-triggered fragmentation 
function as a function of momentum fraction zt in Au-I-Au 
collisions at RHIC, with b = 4.5 fm. 



PHENIX Prehm., p.^.^^ = 3-5 GeV, ly^^, y^l < 0.35 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 




PT,h = 


3-5 GeV 


- - PT,h = 


5-7 GeV 


■ - ■ PT,h = 


7-9 GeV 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIG. 9: (Color online) Photon-triggered Iaa as a function of 
momentum fraction zt in Au-|-Au collisions at RHIC, with 
b = 4.5 fm. 



tons is drastically reduced in absolute magnitude when 
combined with the other components whose momentum 
anisotropy is either zero or positive. Recall that the ther- 
mal photons, not included in this calculation, will kick in 
as the transverse momentum gets lower than vt ~ 3 — 4 
GeV/c, and that their V2 is positive [13, [6^. A com- 
bination of isolation and centrality cuts will be needed 
to clearly isolate the predicted [b^] negative contribution 

So far, we have just considered single jets and single 
photons. Now, results for the correlation between back- 
to-back high-pT trigger photons and associated hadrons 
will be discussed. In Fig. [3 the photon-triggered frag- 
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mentation function Dpp{zT) in P+P collisions at RHIC 
is plotted as a function of momentum fraction for 
three different associated hadron px ranges: = 3- 
5 GeV, 5-7 GeV and 7-9 GeV. The experimental data 
are taken from PHENIX [g^], and we have chosen those 
data points with highest transverse momentum for asso- 
ciated hadrons {pf^ = 3-5 GeV) such that fragmentation 
might be dominant for hadron production. We find that 
the slopes of the photon-triggered fragmentation func- 
tion Dpp(zt) are slightly different for the three hadron 
ranges. We may trace this difference back to the different 
momentum (fraction) dependence of initial parton distri- 
bution functions (PDF) for quarks and gluons (a steeper 
slope for gluon PDF than quark PDF as a function of mo- 
mentum fraction) . As we increase the momenta of initial 
partons (and therefore Bjorkcn x) more quarks appear in 
the initial state, so the relative importance of qq gj 
increases compared to qg qj. As a consequence, the 
photon-triggered fragmentation function Dpp{zT), which 
is an average over quark and gluon fragmentation func- 
tions weighted with their fractions, will become steeper, 
since the gluon fragmentation function is steeper than 
that for quarks. 

In Fig. [HI we show our results for the photon-triggered 
fragmentation function Dj^j^^zt) in mid-ccntral Au-f Au 
collisions (6 = 4.5 fm) as a function of momentum frac- 
tion ZT, compared with 0-20% data from PHENIX [H]. 
The dependence of Daa{zt) on associated hadron pT 
is similar to Dpp{zT) in the low zt regime because the 
associated hadrons are mostly produced from those jets 
triggered by direct photons (see Fig. [T5|) . However, this 
difference tends to decrease at high zt due to the presence 
of jet-medium interaction in Au-|-Au collisions. Since jet- 
plasma photons and fragmentation photons assume more 
significance to photon production at lower (see Fig.S]), 
a larger enhancement is obtained for smaller associated 
hadron pT- 

To further quantify the nuclear medium effect on the 
photon-triggered hadron production, we may take the ra- 
tio of Daa to Dpp and calculate the nuclear modification 
factor Iaa for photon-triggered hadron production. The 
result is shown in Fig. [H] and compared with experimental 
measurement by PHENIX. At low zt, Iaa falls with in- 
creasing Zt due to the dominance of direct photons, but 
it flattens out at higher zt due to the effect of jet-medium 
photons and fragmentation photons (see Fig. [T5|) . 

Also, we compare our calculation with STAR measure- 
ments [g^- In Fig. [ini we show the photon-triggered 
fragmentation function Daa{zt) as a function of mo- 
mentum fraction zt in central Au-|-Au collisions (6 = 
2.4 GcV) at RHIC. The trigger photon momenta have 
been chosen to be pip = 8-16 GeV. The theoretical curve 
of photon-triggered fragmentation function Daa{zt) hi 
central Au-|-Au collisions agrees with experimental data 
quite weU. The data for peripheral (40-80%) Au-|-Au 
collisions are also shown in the figure for the interest of 
comparing with the p-|-p calculation. For completeness, 
the nuclear modification factor Iaa for photon-triggered 
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FIG. 10: (Color online) The photon-triggered fragmentation 
functions as a function of momentum fraction zt in central 
Au-f An collisions and p-f-p collisions (peripheral Au-I-Au col- 
lisions) at RHIC. 
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FIG. 11: (Color online) Photon-triggered Iaa for hadrons as 
a function of hadron momentum pr or momentum fraction 
Zt in Au-|-Au collisions at RHIC, with b = 2.4 fm. 

hadron production is plotted in Fig. [TT] function 
of momentum fraction zt and hadron momentum pi^. 
While Iaa as a function of zt is similar to Fig. [9l Iaa is 
a rather flat function of p^ because flxing p^ with a wide 
range of p^ is actually an average over a variety of zt- 

We also study the medium-size dependence of photon- 
hadron correlations by calculating the photon-triggered 
hadron production in non-central collisions. By integrat- 
ing over the transverse momenta of both trigger photons 
and associated hadrons, we study the centrality depen- 
dence of the per-trigger yield of photon-triggered hadrons 
in Au+Au collisions at RHIC. This result is plotted in 
Fig. 1121 where the experimental measurements arc from 
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FIG. 12: (Color online) The per-trigger yield for photon- 
tagged hadrons in Au+Au collisions at RHIC as a function 
of centrality. The four points in each theoretical curve corre- 
spond to four impact parameters, 2.4 fm , 4.5 fm, 6.3 fm and 
7.5 fm. 



STAR [69| . In the calculation, we have chosen four differ- 
ent impact parameters, 2.4 fm, 4.5 fm, 6.3 fm and 7.5 fm, 
which correspond to four points in each theoretical curve. 
We do not perform the calculation for very peripheral col- 
lisions because the assumption of a thcrmalizcd medium 
essential for a hydrodynamical treatment is no longer ful- 
filled. Again, the trigger photon px has been chosen to 
be = 8-16 GeV, and two different hadron pT ranges 
(pt^ = 4„6 GeV and 6-8 GeV) are studied. We find that 
the centrality dependence of per-trigger yield from the 
calculation is also consistent with current experimental 
measurements. 
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FIG. 13: (Color online) The contributions from different pho- 
ton channels to the photon-triggered fragmentation function 
in p+p collisions and Au+Au collisions at RHIC. 

In the above results, we have taken into account all 



possible sources of high-pr photon production. It would 
be interesting to study how different sources of photons 
contribute the final photon-triggered hadron production. 
To address this issue, we may decompose the per-trigger 
yield of hadrons into different parts, each associated with 
a specific photon source. 



= ^P(p^,srcb5,), 



(17) 



where P(p^, srcjpi^) = P(p^,py, src)/P(py). A similar 
decomposition may be performed for photon-triggered 
fragmentation functions. In Fig. 1131 we show the rela- 
tive contributions from different photon sources to the 
photon-triggered fragmentation function D(zt,P^)) m 
both p+p collisions and central Au+Au collisions {b = 
2.4 fm) if we trigger on a photon on the near side. We 
find that about half of the away-side hadrons at relatively 
small zt are produced from those jets tagged by direct 
photons, while at higher zt, a large amount of away- 
side hadrons come from those jets tagged by jet-medium 
photons and fragmentation photons. Especially, close to 
Zt = 1, where the associated hadrons have almost the 
same amount of transverse momentum as the trigger pho- 
ton, the away side hadron production is dominated by 
those jets tagged by fragmentation photons. Therefore, 
it would be interesting to extend the experimental data 
to higher zt (and to high pt^ such that fragmentation 
dominates hadron production), where jet-medium inter- 
action and fragmentation make a significant contribution 
to photon-hadron correlations. 



VI. CONCLUSIONS 

In this work, we have studied jet energy loss and pho- 
ton production at high pT together in the same frame- 
work, with additional information provided by photon- 
hadron correlations. The energy loss of hard jets travers- 
ing the hot and dense medium is computed by a consis- 
tent incorporation of both induced gluon radiation and 
elastic collisions. The production of high pT photons is 
obtained by taking into account a complete set of photon- 
production channels. Numerical results have been pre- 
sented and compared with experimental measurements 
by employing a fully (3+l)-dimensional hydrodynamical 
evolution model for the description of the thcrmalizcd 
medium created at RHIC. 

Our results illustrate that the magnitude of jet quench- 
ing in relativistic nuclear collisions is sensitive to the in- 
clusion of both radiative and coUisional energy loss. It 
is also found that the interaction between jets and the 
thcrmalizcd medium makes a significant contribution to 
the total photon production at RHIC. Therefore, it is 
important to incorporate all sources of hard photons for 
a full understanding of the correlations between back-to- 
back hard photons and hadrons. In summary, our study 
provides the groundwork to experimentally test our un- 
derstanding of the jet-medium interaction and further 
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extract detailed information about the hot and dense 
medium created at RHIC. 
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